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HIGHLIGHTS 


•  An  extended  equivalent  circuit  battery  model  (EECM)  with  enhanced  accuracy  in  low  state-of-charge  (SOC)  area  is  introduced. 

•  The  EECM  uses  the  surface  SOC  (SOCsurf)  to  reflect  the  solid-phase  diffusion  process. 

•  The  low-SOC-area  voltage  estimation  error  could  be  reduced  by  more  than  50%  through  the  EECM. 

•  An  improved  SOC  estimation  accuracy  in  low-SOC  area  is  resulted  by  the  EECM. 
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In  order  to  predict  the  battery  remaining  discharge  energy  in  electric  vehicles,  an  accurate  onboard 
battery  model  is  needed  for  the  terminal  voltage  and  state  of  charge  (SOC)  estimation  in  the  whole 
SOC  range.  However,  the  commonly-used  equivalent  circuit  model  (ECM)  provides  limited  accuracy  in 
low-SOC  area,  which  hinders  the  full  use  of  battery  remaining  energy.  To  improve  the  low-SOC-area 
performance,  this  paper  presents  an  extended  equivalent  circuit  model  (EECM)  based  on  single¬ 
particle  electrochemical  model.  In  EECM,  the  solid-phase  diffusion  process  is  represented  by  the 
SOC  difference  within  the  electrode  particle,  and  the  terminal  voltage  is  determined  by  the  surface  SOC 
(SOCsurf)  representing  the  lithium  concentration  at  the  particle  surface.  Based  on  a  large-format 
lithium-ion  battery,  the  voltage  estimation  performance  of  ECM  and  EECM  is  compared  in  the  entire 
SOC  range  (0-100%)  under  different  load  profiles,  and  the  genetic  algorithm  is  implemented  in  model 
parameterization.  Results  imply  that  the  EECM  could  reduce  the  voltage  error  by  more  than  50%  in 
low-SOC  area.  The  SOC  estimation  accuracy  is  then  discussed  employing  the  extended  Kalman  filter, 
and  the  EECM  also  exhibits  significant  advantage.  As  a  result,  the  EECM  is  very  potential  for  real-time 
applications  to  enhance  the  voltage  and  SOC  estimation  precision  especially  for  low-SOC  cases. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Due  to  the  excellent  performance  in  driving  economy  and 
environmental  friendliness,  electric  vehicles  (EVs)  have  attracted 
significant  attention,  especially  for  the  pure  electric  vehicles  (PEVs) 
with  zero-emission  capability.  However,  the  PEV  driving  range  is 
not  competitive  compared  with  the  range  of  conventional  vehicles 
owing  to  the  low  energy  density  and  high  cost  of  vehicle  battery, 
which  is  one  of  the  key  drawbacks  from  the  view  of  the  potential  EV 
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consumers  [1].  Additionally,  the  estimation  error  of  PEV  remaining 
driving  range  (RDR)  could  not  be  neglected  owing  to  the  inaccuracy 
of  battery  state  estimation  and  brings  considerable  range  anxiety  to 
the  passengers  [2  .  An  accurate  RDR  depends  on  the  precise  esti¬ 
mation  of  battery  remaining  discharge  energy  (Erde),  which  relies 
on  an  accurate  battery  model  to  calculate  the  terminal  voltage  (Dt) 
and  the  state  of  charge  (SOC).  As  in  Eq.  (1 ),  the  remaining  discharge 
energy  (Erde)  could  be  calculated  by  the  future  average  terminal 
voltage  (Dtiavg)»  the  current  SOC,  and  the  permitted  lower  SOC  limit 
(SOQim,/)  to  prevent  overdischarge,  in  which  /  is  the  future  current 
input  and  Qst  is  the  battery  standard  capacity  [3].  The  SOCiim>/  is 
determined  by  the  onboard  battery  management  strategy  consid¬ 
ering  SOC  estimation  error. 
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Lithium-ion  battery  is  one  of  the  best  candidates  for  PEV  energy 
storage  system  for  its  high  energy  density  and  cycling  lifetime,  and 
a  plenty  of  real-time  modeling  efforts  have  been  made  [4-13]. 
However,  since  the  Li-ion  battery  shows  strong  nonlinearity  in  low- 
SOC  area  (e.g.  at  SOC  lower  than  20%),  the  present  onboard  models 
could  hardly  provide  satisfying  SOC  and  voltage  results  in  low-SOC 
range,  and  an  overestimation  of  SOC  could  lead  to  battery  energy 
exhaustion  before  expectation  and  force  the  PEV  into  limp-home 
status.  As  a  result,  a  certain  part  of  battery  energy  needs  to  be 
reserved  in  real-world  applications  by  technically  setting  a  rela¬ 
tively  large  SOQim,/  (normally  larger  than  10%)  [14,15].  The  battery 
energy  could  hence  not  be  fully  implemented  and  the  PEV  driving 
range  is  further  restricted.  If  the  battery  model  accuracy  in  low-SOC 
area  is  raised,  a  smaller  SOCiim>/  value  could  be  used  and  more 
battery  energy  could  hence  be  utilized  to  obtain  a  larger  driving 
range  without  battery  cost  increase.  To  our  knowledge,  the  low- 
SOC-area  model  accuracy  is  seldom  discussed  in  existing 
publications. 

There  exist  a  number  of  studies  on  lithium-ion  battery  models, 
in  which  the  main  work  lies  in  equivalent  circuit  models  [6-13], 
electrochemical  models  [16-18  ,  and  the  impedance  model  19,20]. 
The  equivalent  circuit  model  (ECM)  reflects  the  battery  perfor¬ 
mance  by  an  electrical  circuit  composed  of  a  set  of  basic  electric 
components  (resistors,  capacitors,  inductors,  etc.),  through  which 
the  voltage  response  under  a  certain  current  input  is  calculated.  The 
model  structure  is  mostly  empirical  based  and  the  parameters  are 
determined  through  battery  discharging/charging  experiments. 
Due  to  its  low  computational  complexity  and  low  difficulty  in 
parameterization  process,  the  ECM  is  frequently  employed  in  on¬ 
board  state  estimation  cases.  The  most  commonly  used  ECM  con¬ 
sists  of  a  pure  resistor  and  several  resistor-capacitor  parallel  circuits 
(RC  components),  such  as  the  one-state  ECM  with  one  RC  compo¬ 
nent  [7,8]  and  the  two-state  ECM  with  two  RC  components  [9-13], 
etc.  Generally,  more  RC  components  in  the  ECM  lead  to  improved 
voltage  and  SOC  estimation  accuracy  at  the  cost  of  larger  calcula¬ 
tion  amount,  and  many  researchers  implemented  the  two-state 
ECM  as  a  balanced  choice  to  validate  the  performance  of  their 
battery  state  estimators.  According  to  the  literature,  the  ECM  could 
provide  relatively  satisfying  voltage  and  SOC  results  in  high  and 
middle  SOC  area,  but  the  model  accuracy  in  low-SOC  area  was 
seldom  mentioned.  The  lowest  SOC  value  was  usually  determined 
not  less  than  10%  in  experimental  validation,  e.g.  40%  by  Hu  et  al. 
[4],  25%  by  Lin  et  al.  [9 , 20%  by  He  et  al.  [6],  and  10%  by  Xiong  et  al. 
[8].  From  our  experience  and  some  of  the  literature  results  [13,19], 
the  precision  of  the  traditional  ECM  usually  becomes  worse  at  SOC 
lower  than  20%  and  decreases  evidently  under  10%  SOC.  As  a  result, 
the  model  inaccuracy  at  low  SOC  hinders  the  ECM  performance  in 
onboard  Erde  estimation,  which  requires  high  model  accuracy  in 
the  whole  SOC  range. 

The  electrochemistry-based  battery  model  is  another  focus  of 
the  present  research,  in  which  terminal  voltage  and  SOC  are 
determined  based  on  the  porous  electrode  structure  and  the  elec¬ 
trochemical  reaction  process,  and  a  better  understanding  of  the 
electrochemical  mechanism  is  provided  compared  with  the 
equivalent  circuit  models.  A  large  number  of  studies  in  this  topic 
are  based  on  the  electrochemical  pseudo-two-dimensional  model 
(P2DM)  developed  by  Doyle  et  al.  [16]  to  describe  the  battery 
macroscopic  behavior  during  operation.  In  Doyle's  model,  the  cell  is 
divided  into  two  dimensions  (i.e.  the  macroscopic  electrode¬ 


separator  dimension  and  the  microscopic  particle  dimension)  and 
the  terminal  voltage  Ut  is  calculated  through  numerical  method  by 
five  partial  differential  equations.  The  large  calculation  load  caused 
by  the  numerical  computation  makes  the  P2DM  hardly  accessible 
in  onboard  cases.  A  series  of  model  simplification  attempts  were 
made  through  different  assumptions  and  approximations  to  reduce 
the  computational  complexity,  and  the  resulted  models  were 
combined  with  real-time  estimation  methods  e.g.  Kalman  filtering 
to  provide  voltage  and  SOC  estimation  results.  By  assuming  the 
current  density  distribution  within  one  electrode  as  uniform  and 
neglecting  the  spatial  difference  in  liquid-phase  potential  and 
liquid-phase  lithium  concentration,  Santhanagopalan  et  al.  [18] 
came  out  with  the  single-particle  model  (SPM)  and  proved  its  ac¬ 
curacy  under  constant-current  discharge  profile.  Based  on  San- 
thanagopalan's  SPM,  Domenico  et  al.  [21]  considered  the  liquid- 
phase  potential  distribution,  Prada  et  al.  [22]  and  Luo  et  al.  [23] 
discussed  the  influence  of  liquid-phase  concentration  distribu¬ 
tion,  and  the  double  layer  effect  was  mentioned  in  Marcicki's  model 
[24].  The  SPMs  present  accurate  results  in  voltage  and  SOC  esti¬ 
mation  in  the  whole  SOC  range,  but  the  increase  of  computational 
load  is  still  not  negligible  compared  with  the  ECM,  and  the  battery 
electrochemical  parameters  could  hardly  be  precisely  determined 
by  automotive  engineers.  Consequently,  the  electrochemical  model 
needs  more  simplification  for  real-world  PEV  applications. 

In  addition,  some  battery  models  are  based  on  the  electro¬ 
chemical  impedance  spectroscopy  (EIS)  method,  through  which  a 
relatively  complex  equivalent  circuit  model  (or  so-called  imped¬ 
ance  model)  is  resulted.  Andre  et  al.  19]  developed  an  impedance 
model  from  EIS  results,  in  which  two  constant  phase  elements 
(CPEs)  and  a  Warburg  component  are  implemented  to  meet  the 
requirements  in  different  current  frequencies.  Li  et  al.  20]  took 
account  of  the  electrolyte  resistance,  the  impedance  of  solid- 
electrolyte  interface  (SEI),  the  double  layer  effect,  and  the  solid- 
phase  diffusion,  and  presented  a  complex  electrochemical  imped¬ 
ance  model  which  perfectly  agreed  with  the  EIS  experimental  re¬ 
sults.  The  impedance  model  provides  better  estimation  results  than 
the  ECM  in  a  wide  frequency  range,  but  the  implementation  of 
complex  electrical  components  (e.g.  the  CPE  and  the  Warburg 
element)  makes  it  impossible  to  transfer  the  model  parameters 
from  the  frequency  domain  into  the  time  domain  directly,  and  the 
impedance  model  is  hence  difficult  to  be  expressed  in  an  iterative 
form  for  real-time  cases.  Similar  to  the  ECM,  the  impedance  model 
also  encounters  problems  in  low-SOC-area  performance  19]. 

The  general  characteristics  of  the  ECM,  the  electrochemical 
model,  and  the  impedance  model  are  summarized  in  Table  1.  As  the 
only  model  type  suitable  for  real-time  cases,  the  ECM  is  chosen  as 
the  basis  of  model  modification  to  improve  the  low-SOC-area 
performance. 

This  paper  aims  at  developing  a  real-time-applicable  battery 
model  with  enhanced  estimation  accuracy  in  low-SOC  area,  through 
which  more  energy  stored  in  the  battery  could  be  used  for  EV  drive. 
Based  on  the  single-particle  electrochemical  model,  this  research 
presents  the  extended  equivalent  circuit  model  (EECM),  in  which 
the  surface  SOC  (SOCsurf)  representing  the  lithium  concentration  at 


Table  1 

Comparison  of  different  types  of  lithium-ion  battery  models. 


Model  structure 

Model  accuracy 
in  low-SOC  area 

Calibration 

effort 

Availability  in 
onboard  cases 

Equivalent  circuit 

Poor 

Low 

Good 

model 

Electrochemical 

Good 

High 

Poor 

model 

Impedance  model 

Poor 

Relatively  low 

Relatively  poor 
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electrode  particle  surface  is  determined  to  describe  the  solid-phase 
diffusion  process.  The  EECM  shows  good  nonlinear  performance 
with  the  help  of  the  SOCsurf  term  and  could  follow  the  voltage 
fluctuation  in  the  entire  SOC  range  very  well.  A  type  of  large-format 
lithium-ion  battery  was  employed  for  experimental  validation,  and 
the  EECM  tends  out  to  provide  largely-improved  performance  in 
low-SOC  area  compared  with  the  commonly-used  ECM.  The 
remainder  of  this  paper  is  structured  as  follows.  Section  2  derives 
the  extended  equivalent  circuit  battery  model  based  on  the 
simplification  of  single-particle  model,  and  the  advantage  of  the 
EECM  in  low-SOC  area  is  theoretically  explained.  Section  3  explains 
the  experimental  setup  and  the  parameter  identification  process. 
The  model  validation  and  comparison  are  discussed  in  detail  in 
Section  4,  and  Section  5  shows  the  real-time  SOC  estimation  results 
through  extended  Kalman  filter.  Section  6  concludes  this  paper. 

2.  Battery  extended  equivalent  circuit  model 

2.1.  Simplified  explanation  of  battery  voltage  response  based  on 
single-particle  model 

From  the  electrochemical  viewpoint,  the  operating  process  of 
lithium-ion  battery  is  explained  as  the  lithium  intercalation/dein¬ 
tercalation  in  electrode  particles  and  the  movement  in  electrolyte 
between  negative  and  positive  electrodes.  The  purpose  of  this 
research  is  to  obtain  a  real-time-realizable  battery  model  through 
simplification  of  the  original  electrochemical  understandings,  and 
the  single-particle  model  (SPM)  is  hence  employed  as  the  basis  of 
the  following  discussion.  The  SPM  treats  the  size  and  dynamics  of 
the  electrode  particles  as  well  as  the  current  density  as  identical 


within  an  electrode,  and  the  negative  and  positive  electrodes  could 
consequently  be  seen  as  two  electrode  particles.  As  in  Eq.  (2),  the 
battery  open  circuit  voltage  (Docv)  is  related  to  the  positive  elec¬ 
trode  potential  Dp(cb,p)  and  negative  electrode  potential  Hn(cb,n), 
which  depends  separately  on  the  bulk  average  lithium  concentra¬ 
tion  of  the  positive  and  negative  electrode  particles  (cb,p  and  Cb,n). 
Note  that  the  variables  related  to  lithium  concentration  in  this  work 
indicate  the  solid-phase  concentration  without  special  statements. 
The  terminal  voltage  Ut  is  more  directly  reflected  by  the  lithium 
concentration  at  the  particle  surface  (csurf>p  and  csurf>n),  i.e.  the 
positive  electrode  surface  potential  Hp(csurfiP)  and  the  negative 
electrode  surface  potential  Hn(csurfin),  which  the  potential  sub¬ 
traction  is  noted  as  ADsurf  in  Eq.  (3).  The  general  relationship  of 
ALlgurf  and  Docv  during  operation  could  be  explained  by  Eq.  (4). 
When  the  battery  is  after  a  long-time  rest,  the  particle  surface 
concentration  equals  the  bulk  average  concentration  (cb,p  and  Cb,n) 
and  ALlsurf  equals  Docv- 

C)CV  =  17p(cb,p)  -  Un(cb  n)  (2) 

^f^surf  =  kfp^Csurfp^  —  k^n^Csurfn^  (3) 

^f^surf  =  k^OCV  "T  (^p(^surf,p^)  —  tfp^Qyp^)  "f  (^n  (cSurf,n) 

-Un(cb,n))  (4) 

The  battery  discharge  process  is  explained  in  Fig.  1.  Here  we 
mainly  consider  the  lithium-ion  battery  with  carbon-based 


Negative  Double  SEI  Electrolyte  and  Double  Positive 
Electrode  Layer  Seperator  Layer  Electrode 

Particle  Particle 


e 


(1) ->(2):  Solid-phase  diffusion  in  negative  electrode  particle 

Voltage  change :  Un(csurf>n)-Un(cb>n) 

(2) ->(3):  Charge  transfer  in  double  layer  (anode) 

Voltage  change :  CDLi/? 

(3) ->(4):  Current  flow  through  SEI  film 

Voltage  change :  CSei 

(4) ->(5):  Current  flow  through  seperator  and  diffusion  in  electrolyte 

Voltage  change :  Co,  CED 

(5) ->(6):  Charge  transfer  in  double  layer  (cathode) 

Voltage  change :  UoL  p 

(6) ->(7):  Solid-phase  diffusion  in  positive  electrode  particle 

Voltage  change:  Up(csurfp)-Un(cb>p) 


Fig.  1.  Explanation  of  lithium-ion  battery  discharge  process  based  on  single-particle  model. 
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negative  electrode,  in  which  the  solid-electrolyte  interface  (SEI) 
film  on  the  anode  surface  needs  consideration.  Lithium  diffuses  in 
anode  particle  to  the  particle  surface  where  the  charge  transfer 
occurs,  moves  through  the  SEI  film  into  electrolyte  [25],  flows  in 
the  liquid  phase  through  the  separator  to  the  cathode  charge- 
transfer  interface,  intercalates  into  the  cathode  particle,  and 
starts  the  cathode  solid-phase  diffusion.  Each  procedure  results  in 
a  corresponded  voltage  change,  which  makes  the  battery  terminal 
voltage  deviating  from  Lfocv-  The  difference  between  surface 
concentration  csurf  and  bulk  concentration  Cb  indicating  the  solid- 
phase  diffusion  results  in  the  gap  between  ADsurf  and  Lfocv 
((Jn(Csurf,n)_lTn(cb,n)  and  Lfp(csurf  p)-Lfp(Cbjp),  noted  later  as  LfsD.n  and 
I/sd.p).  and  the  double  layer  effect  in  anode  and  cathode  leads  to 
the  voltage  changes  ElDL,n  and  DDl,p.  The  voltage  drops  Usei  and  Uo 
are  caused  by  the  current  flow  through  the  SEI  film  and  the 
electrolyte,  and  the  effect  of  contact  resistance  and  current  col¬ 
lector  resistance  is  also  combined  in  Uo  for  convenience.  UEd 
represents  the  voltage  drop  owing  to  liquid-phase  diffusion.  The 
battery  voltage  output  could  then  be  expressed  in  Eq.  (5).  It  is 
noticed  that  the  voltage  change  in  SEI  film  is  a  high-frequency 
term  with  the  time  constant  lower  than  0.1  s  [20  ,  and  the  film 
could  be  seen  as  a  pure  resistor  for  the  low-frequency  EV  driving 
conditions,  in  which  Usei  is  merged  with  Uo.  The  voltage  changes 
on  the  double  layer  capacitance  (UDLtn  and  UDEiP)  are  later  analyzed 
as  transient  processes.  The  voltage  drop  caused  by  electrolyte 
diffusion  needs  to  be  considered  only  under  very  high  current  rate 
and  could  only  be  well  calculated  by  more  than  seven  parameters 
[23],  which  largely  raise  the  model  complexity.  The  term  UEd  is 
consequently  here  neglected  and  the  terminal  voltage  could  be 
written  as  Eq.  (6). 

Ut  -  Lfocv  -  UsD,n  -  ^SD,p  ~  Uo  ~  ^SEI  —  ^DL,n  ~  ^DL,p  —  ^ED 

(5) 

=  Lfocv  -  fisD,n  -  ^SD,p  —  Uq  —  UDE  n  -  UDE  p  (6) 

According  to  the  SPM,  the  battery  state  of  charge  (SOC)  is  rep¬ 
resented  by  the  average  bulk  concentration  in  the  electrode  particle 
(Cb).  For  the  negative  electrode,  the  bulk  concentration  when  the 
battery  is  full  (SOC  equals  100%)  is  noted  as  Cb,n,ioo%.  and  the  con¬ 
centration  for  an  empty  battery  (SOC  equals  0%)  is  written  as  Cb,n,o%- 
As  a  result,  the  present  SOC  could  be  calculated  by  Cb,n  or  the 
positive  electrode  concentration  CbiP  as  in  Eq.  (7),  and  the  open 
circuit  voltage  Lfocv  is  hence  related  to  the  macroscopic  SOC. 

gQC  __  Cb,n  ~  cb,n,0%  _  cb,p  ~  cb,p,0%  ^ 

cb,n,100%  “  cb,n,0%  cb,p,100%  ~~  cb,p,0% 

Known  from  Eqs.  (4)  and  (6),  Ut  and  ADsurf  are  more  directly 
related  to  the  surface  concentration  csurf  rather  than  the  bulk 
concentration  Cb  indicating  macroscopic  SOC,  and  the  imple¬ 
mentation  of  SOC  could  not  reflect  the  battery  dynamics  straightly. 
Here  we  define  the  term  surface  state  of  charge  (SOCsurf)  to  reflect 
the  particle  surface  lithium  concentration  in  Eq.  (8).  The  surface 
SOCs  of  anode  and  cathode  are  in  direct  relation  to  the  present 
potentials  at  the  surface  of  positive  and  negative  electrodes,  and  the 
terminal  voltage  could  consequently  be  represented  by  SOCsurf  in 
Eq.  (9).  For  the  battery  management  system  (BMS)  developers 
without  specific  electrochemical  techniques,  separate  calibration  of 
the  negative  and  positive  electrode  potentials  is  very  inconvenient. 
To  make  our  presented  model  more  suitable  for  real-world  appli¬ 
cations,  a  single  SOCSUrf-potential  curve  is  used  to  approximate  the 
subtraction  of  positive  surface  potential  Hp(SOCsurfjP)  and  negative 
surface  potential  Hn(SOCsurf>n),  i.e.  AHsurf.  The  SOCsurf-AHsurf  rela¬ 
tionship  could  actually  be  represented  by  the  SO C-Hqcv  curve  with 


SOCsurf  as  the  input,  and  Eq.  (9)  is  hence  rewritten  as  Eq.  (10)  while 
AHsurf  is  later  represented  by  Hocv(SOCsurf)  as  in  Eq.  (11).  The 
SOC-Hocv  calibration  is  very  common  in  present  BMS  development 
and  the  term  SOCsurf  could  then  be  applied  onboard  without 
difficulty. 

SOCsurfn  =^ILJ1-S0C.  SOCsurfiP=^-SOC  (8) 

cb,n  cb,p 

—  (Up  ^SOCsurf  —  Un  (SOCsurf —  Uo  —  hfp)L,n  —  ^DL,p 

(9) 

Ut  =  Uocv(sOCsurf)  -Uo-  UDUn  -  t/DLp  (10) 

A^surf  =  Up  (SOCsurf  p)  —  (in  (SOCsurf  n)  =  L/qcv  (SOCsurf) 

(11) 

During  battery  operation,  the  macroscopic  SOC  is  achievable  by 
current  integration  method  while  the  surface  SOC  could  not  be 
directly  determined.  As  a  result,  the  difference  between  SOCsurf  and 
SOC  (related  to  the  difference  of  csurf  and  Cb  in  solid-phase  diffusion 
process  following  the  Fields  second  law)  requires  discussion  to 
determine  the  terminal  voltage.  Considering  the  double  layer  effect 
[26],  the  current  density  j  in  the  electrode  is  divided  into  the  solid- 
phase-diffusional  current  density  jSD  effective  for  the  solid-phase 
diffusion  and  the  double-layer  current  density  jDL,  and  the  three 
current  densities  (J,  jSd  and  jDL)  correspond  separately  to  the 
macroscopic  current  J,  solid-phase-diffusional  current  Isd  and 
double  layer  current  Jdl  as  Eq.  (12).  The  solid-phase-diffusional 
current  Isd  is  seen  for  simplicity  as  in  one-state  relationship  with 
the  current  /  with  the  time  constant  rDL,  as  in  Eq.  (13),  and  the 
voltage  drop  in  the  double  layer  (Ddl)  is  resulted  from  the  resis¬ 
tance  Rqt  indicating  the  resistance  effect  of  double  layer  in  Eq.  (14). 
In  steady  state  under  constant  current  input,  the  difference  of  csurf 
and  Cb  (noted  as  Acstat,cc)  is  expressed  in  Eq.  (15)  [27],  in  which  the 
concentration  difference  parameter  /<sd,c  is  defined  with  the  unit 
mol  •  m-3  A-1.  The  particle  radius  Rs,  active  material  volume  fraction 
£“s,  electrode  thickness  <5,  electrode  plate  area  A,  and  solid-phase 
diffusion  coefficient  Ds  are  all  constant  for  a  specific  battery  un¬ 
der  constant  temperature,  so  /<sd,c  is  seen  as  constant  and  Acstat,cc  is 
proportional  to  current  input.  The  concentration  difference  Ac 
changes  with  time  in  dynamic  current  cases  and  is  calculated  by  the 
finite  difference  method  in  Ref.  [21  ,  and  the  Ac  variation  under  a 
current  pulse  was  approximated  as  a  first-order  process  for  each 
electrode  in  Ref.  [28  ,  as  in  Eq.  (16)  with  tsd  the  time  constant.  It  is 
to  be  noted  that  the  concentration  difference  Ac  is  related  to  a 
single  electrode,  and  the  difference  between  SOC  and  SOCsurf 
(ASOC)  reflects  the  influence  of  both  electrodes.  The  ASOC  is  hence 
determined  by  two  first-order  processes  in  Eq.  (17),  in  which  /<sd,i 
and  /<sd,2  are  the  SOC-difference  parameters  and  tsd.i  and  tsd,2  are 
the  time  constants. 


j  =Jdl  +7sd^  =  hi  +  *sd  (12) 

Isd  =  J*  [1  -  exP(  -  T/ tDl)]  (13) 

Udl  =  Isd^ct  =  JKcr[l  ~  exP(  ~  T/t"dl)]  (14) 


ACstat,cc  —  Q)  ^surf 


^slSD  _  ^s^SD  _  i  j 

15 DsesF  15 DsesFA5  SD  c  SD 


(15) 


M.  Ouyang  et  at  /  Journal  of  Power  Sources  270  (2014)  221—237 


225 


AC  =  /<SD,  JsD  •  [1  -  exP(  -  f/r SD)] 


(16) 


Upo  =  ^SD,ECM  +  ^DL,ECM 


(22) 


ASOC 


:  ^SD,l^SD,l  '  [1  - 
+  ^SD,2^SD,2  '  [1 


exp(  -  t/rSD,l)' 
-  exp(  -  t/rsD,2) 


(17) 


2.2.  Battery  extended  equivalent  circuit  model  and  model 
comparison 

In  this  section,  the  voltage  calculation  process  of  the  commonly- 
used  battery  models  is  evaluated  according  to  the  SPM-based 
analysis  in  Section  2.1,  and  the  extended  equivalent  circuit  model 
(EECM)  is  introduced  with  the  enhanced  estimation  accuracy 
compared  with  the  traditional  ECM.  The  onboard-frequently- 
employed  ECM  consists  of  a  pure  resistor  and  several  RC  compo¬ 
nents,  and  the  following  discussion  deals  with  the  two-state  ECM 
as  shown  in  Fig.  2,  in  which  the  terminal  voltage  is  expressed  as  Eqs. 
(18)  and  (19).  Referencing  the  electrochemistry-based  explanation 
in  Eq.  (6),  the  pure  resistor  R0  in  the  ECM  reflects  the  instantaneous 
voltage  drop  Do.  The  voltage  source  is  looked  up  by  the  macroscopic 
SOC  with  the  resulted  Docv(SOC)  representing  the  open  circuit 
voltage  instead  of  the  present  potential  D0cv(SOCsurf),  indicating 
that  the  difference  between  Docv(SOC)  and  Docv(SOCsurf)  referring 
to  the  solid-phase  diffusion  (noted  as  Usd  in  Eq.  (20))  needs  specific 
consideration.  The  voltage  drops  owing  to  double  layer  effect  (DoL,n 
and  1/dl.p)  could  be  determined  by  two  one-state  processes  refer¬ 
ring  to  Eq.  (14)  and  hence  included  in  the  polarization  voltages  L7po,i 
and  Upo  2,  noted  as  DDLiniECM  and  Ddl>PiECm  respectively.  As  the  only 
transient  components  in  the  ECM,  the  two  RC  components  also 
reflect  the  solid-phase-diffusional  voltage  drop  Usd  (the  term  re¬ 
flected  in  the  ECM  noted  as  Dsd.ecm)  to  a  certain  extent,  as  in  Eq. 
(21).  In  Eq.  (22)  for  convenience,  the  polarization  voltage  terms 
(Lfpo.i  and  t/po,2 )  are  considered  together  as  Upo,  as  well  as  the 
double-layer  voltage  drops  ElDL,n,ECM  and  DDLiPiECM  (considered 
together  as  DDl,ecm). 


Ut  =  Uocv(SOC)  -  I  Ro  -  UpoA  -  Upo, 2 


(18) 


^po,i  —  ^po,i '  exp(  t/Tpo,i)_  5  1  —  1,2 


(19) 


Usd  =  Lfocv  ( SOC )  -  fiocv  ( SOCsurf 


(20) 


^OCV,ECM  ( SOCsurf )  =  Lf0cv  SOC  +  Usd  ,ECM 


(21) 


Compared  with  the  electrochemistry-based  SPM  discussed  in 
Section  2.1  shown  in  Fig.  3(a),  the  ECM  calculates  the  solid-phase- 
diffusional  voltage  change  USd  and  double-layer-based  voltage 
changes  Udl  in  a  much  simpler  way,  as  illustrated  in  Fig.  3(b).  For 
the  simplified  electrochemical  process  in  Fig.  3(a),  the  solid-phase- 
diffusional  current  Isd  and  the  UD i  are  in  one-state  relationship 
with  the  input  current  7  according  to  Eqs.  (13)  and  (14).  The  dif¬ 
ference  between  SOC  and  SOCsurf  (ASOC)  related  to  the  solid-phase 
diffusion  is  also  one-state  related  to  7Sd-  To  determine  Us d,  the  open 
circuit  voltage  Docv(SOC)  and  the  Docv(SOCSUrf)  are  separately 
looked  up  by  SOC  and  SOCsurf  from  the  SOC-Docv  relationship.  Due 
to  the  strong  nonlinearity  of  the  SOC-Docv  curve  especially  in  low- 
SOC  area,  Usd  could  not  be  proportionally  related  to  ASOC.  As  a 
result,  the  solid-phase-diffusional  Usd  is  determined  through  a 
second-order  process  and  a  nonlinear  mapping  relation.  In 
contrast,  the  solid-phase-diffusional  voltage  change  Dsd.ecm  for  the 
ECM  is  directly  in  one-state  relationship  with  the  macroscopic 
current  as  part  of  the  polarization  voltage  Dpo.  As  a  result,  the  ECM- 
determined  Dsd.ecm  may  be  inaccurate  due  to  the  following  rea¬ 
sons.  Firstly,  the  first-order  process  in  ECM  is  used  to  simulate  the 
second-order  relationship  in  SPM,  which  lowers  the  voltage  accu¬ 
racy.  Secondly,  Dsd.ecm  and  Ddl,ecm  are  expressed  together  by  Dpo. 
Since  the  time  constants  of  the  double  layer  and  the  solid-phase 
diffusion  processes  are  different,  a  compromise  must  be  made  to 
obtain  the  rpo.  And  the  Dsd.ecm  could  not  reflect  the  nonlinear 
relationship  of  Usd  and  ASOC,  which  might  be  the  key  drawback  of 
the  traditional  ECM. 

In  order  to  explain  the  necessity  to  consider  the  SOC-Docv 
nonlinearity,  Fig.  4  discusses  the  determination  of  Docv(SOCsurf) 
under  different  current  rates  as  an  example.  As  a  first  order  process 
in  Eq.  (19),  the  polarization  voltage  Dpo  at  time  point  t  is  enlarged  or 
reduced  by  a  factor  number  a  when  the  input  current  7  is  multiplied 
by  a,  i.e.  Upo(aI,  t)  =  aDpo(7,  t),  and  the  ECM-based  Dsd.ecm  is 
assumed  to  follow  the  same  pattern  Usd,ecm{ciI ,  t)  =  oDsd.ecmU,  t). 
The  similar  relationship  is  also  for  the  second-order  ASOC  in  SPM, 
i.e.  ASOC(a/,  t)  =  aASOC(7,  t).  In  Fig.  4,  the  open  circuit  voltage  is 
determined  by  the  macroscopic  SOC  as  Docv(SOC),  and  the  present 
ADsurf  follows  the  SO C-D0Cv  curve.  At  a  certain  time  point  t,  the 
surface  SOC  under  a  relatively  low  current  7/  is  noted  as  SOCsurf/, 
with  the  potential  as  Docv(SOCsurf,/).  The  difference  between  SOC- 
surfj  and  macroscopic  SOC  is  ASOC/.  Assuming  a  perfect  identifica¬ 
tion  result  for  the  ECM,  in  which  the  solid-phase-diffusional  voltage 


(a)  Simplification  of  SPM 


soc- 


/ 


1  ^DLS 


hn !  ^sd 


!  1 


SOC 


SOC-Cocv 
relationship 


£/ocv(SOC) 
Usd 


surf 

► 


ASOC 


- -Spc - 'Cocv(SOCw/) 


-W  R 


CT 


Cdl 

- ► 


(b)  ECM 


I 

R,, 

Upo—  U dl,ecm+  U sd,ecm 

- w* 

1 - 

+ 

O 

Fig.  3.  Calculation  of  the  solid-phase  diffusional  voltage  change  Usd  and  the  double- 
layer-based  voltage  changes  t/DL  based  on  different  models:  (a)  based  on  simplified 
explanation  of  single-particle  model;  (b)  based  on  equivalent  circuit  model. 
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drop  Usd, ecm, i  equals  the  subtraction  of  the  open  circuit  voltage 
Uocv(SOC)  and  the  Hocv(SOCSUrf,/)-  According  to  Eqs.  (20)  and  (21), 
the  ECM-estimated  H0cv,ECM(SOCsurfj/)  equals  El0cv(SOCsurfi/).  How¬ 
ever,  for  another  current  rate  7h  with  7h  =  bli  ( b  is  the  factor),  the 
estimated  result  under  7/  would  lead  to  a  certain  error.  From  the 
above  analysis,  ASOCh  (the  difference  between  SOCSUrf,h  and 
macroscopic  SOC)  at  time  point  t  equals  bASOC/.  For  the  ECM,  the 
diffusional  voltage  drop  Hso.ECM.h  equals  bHsD,ECM,h-  The  ECM- 
estimated  Hocv,ECM(SOCsurfih)  could  hence  be  determined  geomet¬ 
rically  as  the  intersection  of  the  extrapolation  of  the  blue-dashed 
AB  and  the  y-direction  line  of  SOCsurf  h,  i.e.  point  C,  which  results 
in  a  large  deviation  from  the  real  Hocv(SOCsurf  h)  (point  D).  For  an 
arbitrary  current  input,  the  ECM-estimated  value  is  on  the  extrap¬ 
olation  or  interpolation  of  the  AB  line,  meaning  that  the  ECM  could 
only  provide  linear  approximation  of  the  voltage  drop  in  solid- 
phase  diffusion.  When  the  SOC-Hocv  curve  shows  good  linearity, 
the  ECM-determined  value  deviates  not  much  from  the  real 
Hocv(SOCsurf).  When  the  SOC-Hocv  curve  is  very  nonlinear  (e.g.  the 
case  in  Fig.  4),  the  gap  between  H0cv,ECM(SOCsurf)  and  the  real  po¬ 
tential  is  very  noteworthy. 

For  lithium-ion  batteries,  good  linearity  is  generally  shown  for 
the  SOC-Hocv  curve  in  high  and  middle  SOC  range,  and  conse¬ 
quently  the  ECM-determined  voltage  illustrates  high  accuracy.  In 
low-SOC  range,  the  open  circuit  voltage  drops  very  rapidly  and 
nonlinearly  as  SOC  decreases,  leading  to  a  large  voltage  estimation 
error.  For  this  reason,  the  ECM  could  not  provide  satisfactory  esti¬ 
mation  results  in  low-SOC  cases.  It  is  also  the  case  for  the  imped¬ 
ance  model  with  more  complex  structure.  Although  the  description 
of  the  solid-phase  diffusion  is  more  accurate  indicating  preciser 
SOCsurf  values,  there  is  still  no  reflection  of  the  SOC-Hocv  curve 
nonlinearity  between  the  voltage  drop  HSd,im  and  the  ASOC.  For  the 
electrochemical  models,  e.g.  the  P2DM  and  the  SPM,  since  the 
solid-phase  diffusion  process  is  specifically  concerned  and  the  cell 
potential  is  determined  by  the  surface  lithium  concentration 
looking  up  the  electrode  potential  relationship,  the  problem  of 
linear  approximation  is  avoided,  and  the  low-SOC  accuracy  is  also 
guaranteed  [18].  But  the  electrochemical  models  depend  on  too 
many  electrochemical  parameters,  which  could  hardly  be  deter¬ 
mined  or  identified  in  real-world  cases.  And  the  high  computa¬ 
tional  complexity  also  makes  these  models  not  applicable  in 
onboard  BMS  chips. 

The  purpose  of  this  paper  is  to  present  an  accurate  onboard- 
applicable  battery  model  which  represents  the  electrochemical 


process  better  than  the  traditional  ECM.  As  discussed  above,  the 
ECM  could  reflect  many  of  the  electrochemical  terms  in  Eq.  (6)  very 
well,  including  the  instantaneous  voltage  change  Ho  by  the  resistor 
R0  and  the  double-layer-based  voltage  changes  (HDL,n  and  HDl,p)  by 
two  RC  components.  But  the  Hsd  arising  from  the  solid-phase 
diffusion  could  only  be  linearly  approximated  in  the  ECM  as  in 
Fig.  4,  which  brings  unignorable  errors  in  the  nonlinear  period  of 
the  SOC-Hocv  relationship.  As  a  result,  this  research  seeks  to 
introduce  a  component  specifically  reflecting  the  solid-phase 
diffusion  instead  of  the  RC  elements  used  in  the  traditional  ECM. 
The  surface  SOC  is  firstly  estimated  through  the  macroscopic  SOC 
and  the  ASOC,  and  the  present  AHsurf  is  directly  looked  up  as 
Hocv( SOCsurf).  The  introduction  of  SOCsurf  in  the  modified  model 
avoids  the  linear  approximation  through  the  open  circuit  voltage 
Hocv(SOC)  and  the  solid-phase-diffusion  voltage  drop  Hsd,  and  in 
that  case  the  shape  of  the  nonlinear  SOC-Hocv  curve  could  always 
be  followed.  According  to  Eq.  (17),  the  term  ASOC  is  achievable  by 
two  first-order  processes,  and  could  be  expressed  by  two  ASOC 
components  (later  noted  as  DSOC  component).  The  modified  model 
is  named  as  the  extended  equivalent  circuit  model  (EECM)  with  the 
schematic  shown  in  Fig.  5.  The  EECM  consists  of  a  pure  resistor,  two 
RC  elements,  two  DSOC  components,  and  a  voltage  source.  The  RC 
components  correspond  respectively  to  the  double-layer  voltage 
change  in  negative  and  positive  electrode  (noted  as  HDl,i  and  HDl,2, 
as  in  Eq.  (23)),  with  the  solid-phase-diffusional  current  7sd,i  and 
7sd,2  calculated  in  Eq.  (24).  The  DSOC  components  are  also  one-state 
but  with  the  solid-phase-diffusional  current  terms  as  input  and  SOC 
offset  as  output  in  Eq.  (25).  The  output  ASOCi  and  ASOC2  is  then 
combined  with  the  macroscopic  SOC  to  determine  the  SOCsurf  in  Eq. 
(26),  which  is  then  employed  by  the  voltage  source  to  obtain  the 
Hocv(SOCsurf)  directly.  Consequently,  the  terminal  voltage  of  EECM 
is  determined  as  in  Eq.  (27).  Referencing  the  electrochemical  terms 
in  Eq.  (6),  Table  2  compares  the  EECM  and  the  traditional  ECM. 


exp(-t/rDLl)], 


(23) 


1  =  1,2 


(24) 


ASOC, 


ksDikoj 


exp(-  t/rSEM)], 


i  =  1,2 


(25) 


SOCsurf  =  SOC  -  ASOC!  -  ASOC2  (26) 

Ut  =  Uqcv  (sOCsurf )  -IR0-  UDU  -  UDL  2  (27 ) 
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Table  2 

Comparison  of  ECM  and  EECM  in  terminal  voltage  determination. 


Electrochemical  terms 

Related  terms  in  ECM 

Related  terms 

in  EECM 

Instantaneous  voltage 

Resistor  R0 

Resistor  R0 

change  ( U0 ) 

Double-layer-based 

RC  components 

RC  components 

voltage  change  (UDL) 

Solid-phase-diffusional 
voltage  change  ( USD ) 

RC  components 
(linear  approximation) 

DSOC  components 

Voltage  source  output  (Uqcv) 

Macroscopic  SOC 

Surface  SOC  (SOCsurf) 

Compared  with  the  traditional  ECM,  only  two  one-state  DSOC 
elements  are  added  in  the  EECM,  resulting  in  relatively  limited 
increase  in  model  complexity  and  a  promising  implementation  in 
real-world  BMS  applications.  The  introduction  of  SOCsurf  indicates 
the  direct  determination  of  ADsurf  in  dynamic  operations  and  the 
model  could  hence  reflect  the  nonlinearity  of  potential  variation.  In 
that  case,  an  evident  improvement  in  low-SOC-area  performance 
could  be  made  by  the  EECM  in  contrast  to  the  ECM. 

3.  Experimental  and  model  parameter  identification 

3.1.  Experimental  setup 

In  this  research,  a  type  of  commercial  prismatic  lithium  nick¬ 
el-cobalt-manganese  oxide/graphite  (NCM-G)  battery  with  20-Ah 
nominal  capacity  was  tested,  in  which  the  NCM-G  system  acts  as  a 
good  candidate  in  the  PEV  application  for  its  high  specific  energy 
and  good  cycle  life  performance  [29].  The  basic  battery  parameters 
are  listed  in  Table  3.  Note  that  the  C-rate  is  related  to  the  nominal 
capacity  Qnom  (20  Ah),  and  the  SOC  value  is  calculated  according  to 
the  standard  capacity  Qst-  The  charge  and  discharge  process  was 
executed  on  a  Neware  BTS-3000  test  bench  with  current 
range  -100  to  100  A.  A  DGBELL  BE-HL-150M3  temperature 
chamber  (temperature  range  -40  to  150  °C,  temperature  error 
within  ±0.3  °C)  was  used  to  limit  the  temperature  fluctuation.  The 
chamber  was  set  at  25  °C  in  this  research,  and  the  actual  temper¬ 
ature  variation  in  all  the  experiments  was  very  limited  (within 
24.8  °C  and  26.2  °C),  indicating  the  temperature  influence  on  bat¬ 
tery  parameters  as  negligible.  For  real-world  applications  with 
large  temperature  range,  a  parameter  calibration  process  under 
different  temperatures  as  well  as  an  electro-thermal  model  is 
needed,  which  is  beyond  the  main  scope  of  this  work  and  will  be 
discussed  in  the  future. 


The  model  parameterization  is  based  on  a  set  of  experiments 
modified  from  the  hybrid  pulse  power  characterization  (HPPC)  test 
[30].  The  original  HPPC  test  with  10  s  discharge,  40  s  pause  and  10  s 
charge  was  designed  for  the  hybrid  electric  vehicle  application,  and 
the  pulse  profile  is  modified  here  by  lengthening  the  discharge 
period  to  30  s  to  fit  the  PEV  case,  as  the  current  profile  shown  in 
Fig.  6  with  discharge  current  40  A  (2  C)  and  charge  current  30  A 
(1.5  C).  The  constant-current  discharge  period  (0.5  C  discharge  for 
280  s,  for  a  total  of  5%  SOC  change  in  one  current  profile)  after  the 
charge  pulse  is  also  included  to  guarantee  the  resulted  parameters 
suitable  for  the  profiles  with  large  current  variation  as  well  as  those 
with  stable  current.  With  a  5%  SOC  gap,  the  test  profile  was 
executed  21  times  from  100%  to  0%  SOC  to  examine  the  SOC  in¬ 
fluence  on  parameter  identification  results,  and  implemented  on 
two  cells  (namely  cell  1  and  cell  2)  to  check  the  model  performance 
considering  cell  inconsistency.  The  standard  capacity  Qst  is 
approximately  20.8  Ah  for  cell  1  and  20.6  Ah  for  cell  2  during  our 
test. 

In  order  to  obtain  a  detailed  SOC-Docv  relationship  for  the 
EECM  rather  than  interpolating  the  limited  HPPC  test  points,  the 
battery  was  discharged  and  charged  under  0.02  C  (0.4  A),  and  the 
averaged  voltage  from  the  low-rate  discharging  and  charging  pro¬ 
cess  is  treated  as  the  SOC-Hocv  curve,  as  shown  in  Fig.  7,  which 
coincides  well  with  the  HPPC-resulted  points.  It  is  noted  that  the 
curve  in  high  and  middle  SOC  range  is  relatively  linear  and  the 
period  under  15%  SOC  shows  obvious  nonlinearity  and  large  slope. 
This  curve  is  later  implemented  in  the  EECM  and  the  ECM,  in  which 
the  EECM  uses  surface  SOC  to  look  up  the  potential  while  the 
traditional  ECM  uses  SOC  instead. 

The  model  estimation  accuracy  is  validated  depending  on  the 
experiments  under  the  constant-current  discharge  pause  profile 
(DPP)  and  the  dynamic  stress  test  (DST)  profile,  which  would  be 
described  detailedly  in  Section  4. 

3.2.  Parameter  identification  and  results  comparison 

In  this  section,  the  measured  voltage  response  under  HPPC 
profile  is  employed  to  identify  the  parameters  of  the  EECM  under 
different  SOCs,  and  the  parameterization  process  is  executed  on  the 
traditional  ECM  (with  the  two-state  equivalent  circuit  model  as  a 
representative)  for  comparison.  There  exist  a  large  number  of  pa¬ 
rameters  to  be  identified  (six  parameters  for  ECM  and  ten  for  EECM, 
listed  in  Tables  4  and  5),  and  a  robust  and  intelligent  identification 
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Table  3 

Basic  parameters  of  test  battery  according  to  product  specifications. 


40 


Parameter 

Value 

< 

Material 

Li(NixMnyCoi_x_y)02)/graphite 

c 

Nominal  capacity  (Ah) 

20 

2 

Maximum  continuous  discharge 

4  C  (80  A) 

L. 

3 

current  (A) 

O 

Maximum  continuous  charge 

1  C  (20  A) 

current  (A) 

Allowed  voltage  range  (V) 

2.75  to  4.2 

Allowed  temperature  range  for 

-20  to  55 

discharge  (°C) 

Allowed  temperature  range  for 

0  to  45 

charge  (°C) 

Recommended  charging  method 

Constant  current  charge 
(0.5  C,  10  A),  then  constant 
voltage  charge  until  0.05  C  (1  A) 

0 


-20- 


-40 


0 


Step  1 :  40  A  discharge  (30  s) 
Step  2:  pause  (40  s) 

Step  3:  30  A  charge  (10  s) 
Step  4: 10  A  discharge  (280  s) 


50  100 


150  200  250 

Time  (s) 


300  350 


Fig.  6.  Current  profile  for  parameter  identification  of  ECM  and  EECM. 
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Fig.  7.  Comparison  of  calculated  SO C—Uocv  relationship  through  0.02  C  discharge/ 
charge  data  and  measured  SOC -U0Cv  from  HPPC  test. 


method  is  needed  to  obtain  the  optimum  results.  Therefore,  the 
genetic  algorithm  (GA)  developed  by  Holland  [31  is  a  proper 
choice.  The  GA  is  a  global  optimization  method  which  could  find 
the  optimum  for  complex  objective  functions  effectively  and  is 
hence  employed  here  for  the  model  parameterization.  The  model 
voltage  output  follows  Eq.  (27)  for  EECM  and  Eq.  (18)  for  ECM,  and 
the  GA  is  applied  to  determine  the  optimum  parameters  by  mini¬ 
mizing  the  root  mean  square  error  between  voltage  measurement 
and  model-calculated  voltage. 

The  to-be-identified  parameters  of  ECM  and  EECM  are  listed  in 
Tables  4  and  5  referencing  Figs.  2  and  5.  The  ohmic  resistance  values 
in  discharging  and  charging  process  are  separately  determined  and 
noted  as  R0,d  and  R0t c,  and  the  upper  and  lower  bounds  of  the  pa¬ 
rameters  in  GA  calculation  are  empirically  decided. 

The  parameterization  is  based  on  the  experimental  data  of  cell  1, 
and  the  accuracy  of  the  obtained  model  is  examined  by  the  data  of 


Table  4 

Parameters  to  be  identified  in  ECM. 


Parameter 

Description 

Unit 

Ro.d 

Discharging  ohmic  resistance 

Q 

R0,c 

Charging  ohmic  resistance 

Q 

ftpo.l 

Resistance  of  the  first  polarization  process 

Q 

rpo,l 

Time  constant  of  the  first  polarization  process 

s 

^po,2 

Resistance  of  the  second  polarization  process 

Q 

rpo,2 

Time  constant  of  the  second  polarization  process 

s 

Table  5 

Parameters  to  be  identified  in  EECM. 


Parameter  Description  Unit 

R0td  Discharging  ohmic  resistance  Q 

R0'C  Charging  ohmic  resistance  Q 

Rct.i  Resistance  of  the  first  double-layer  process  Q 

rDLj i  Time  constant  of  the  first  double-layer  process  s 

Rct.2  Resistance  of  the  second  double-layer  process  Q 

tdl.2  Time  constant  of  the  second  double-layer  process  s 

/<sd.i  SOC-difference  parameter  of  the  first  ASOC  process  — 

t sd,i  Time  constant  of  the  first  ASOC  process  s 

/<:sd.2  SOC-difference  parameter  of  the  second  ASOC  process  - 

r sd,2  Time  constant  of  the  second  ASOC  process  s 


cell  1  as  well  as  cell  2.  The  parameterization  data  source  and  the 
validation  profile  are  identical  for  the  cell  1  case,  which  is  actually 
the  fidelity  examination  of  the  identification  process.  The  root 
mean  square  errors  (RMSEs)  of  the  model-calculated  at  different 
SOCs  are  shown  in  Fig.  8  as  blue  markers  (in  the  web  version)  for 
ECM  and  red  markers  for  EECM.  Both  models  demonstrate  good 
accuracy  at  SOCs  larger  than  20%,  in  which  the  average  RMSE 
resulted  from  ECM  is  about  2  mV,  while  the  case  for  the  EECM  is 
approximately  1  mV.  In  the  lower-than-20%  SOC  range,  the  differ¬ 
ence  between  the  two  models  becomes  more  obvious  and  the 
EECM  shows  evidently  better  accuracy.  For  instance,  the  error  of 
EECM  at  5%  SOC  is  lower  than  10  mV  while  in  contrast  the  ECM- 
based  RMSE  exceeds  30  mV.  It  is  therefore  illustrated  that  the 
EECM  occupies  a  much  better  performance  than  the  ECM  in  low- 
SOC  area. 

The  model-calculated  voltage  is  compared  with  the  measure¬ 
ment  result  in  Fig.  9.  The  results  at  two  SOC  points  representing  the 
middle-SOC  area  (50%  SOC,  in  Fig.  9a)  and  the  low-SOC  area  (5% 
SOC,  in  Fig.  9b)  are  shown  to  compare  the  fitting  results  of  the  ECM 
and  the  EECM.  At  50%  SOC,  the  voltage  results  from  both  models 
follow  the  experimental  data  tightly,  while  at  5%  SOC  the  EECM- 
determined  voltage  traces  the  nonlinear  variation  of  the  mea¬ 
surement  very  well  to  bring  a  much  smaller  error  than  the  tradi¬ 
tional  ECM.  It  is  consequently  expressed  that  the  EECM  with  SOCsurf 
determination  has  a  certain  advantage  over  the  ECM  in  case  of 
nonlinear  battery  properties  at  low  SOCs. 

It  is  essential  to  check  the  performance  of  identified  model 
parameters  on  different  cells  in  consideration  of  the  cell  inconsis¬ 
tency  [32].  The  model  parameters  identified  from  cell  1  are 
employed  on  the  profile  in  cell  2,  and  the  resulted  RMSEs  are  shown 
in  Fig.  10.  It  is  noticed  that  the  RMSE  values  increases  a  little  bit 
compared  with  Fig.  8  due  to  cell  variance,  while  the  errors  at 
relatively  large  SOC  are  still  very  limited  (3  mV  for  ECM  and  2  mV 
for  EECM).  When  the  SOC  is  under  20%,  high  benefits  could  also  be 
made  by  the  EECM  to  reduce  voltage  error. 

Fig.  11  compares  the  model-determined  voltage  output  and  the 
measured  voltage  on  cell  2,  with  Fig.  11a  at  50%  SOC  and  Fig.  lib  at 
5%  SOC.  The  estimation  results  are  satisfying  for  both  models  at 
50%  SOC  with  only  slightly  RMSE  increase,  and  the  EECM  could 
still  guarantee  its  extraordinary  accuracy  at  low  SOC  on  a  different 
cell.  As  a  result,  the  EECM  is  applicable  for  the  battery  pack  with 
cell  differences  to  enhance  the  low-SOC-area  estimation 
performance. 
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Fig.  8.  Voltage  RMSE  comparison  of  ECM  and  EECM  under  the  parameterization  profile 
at  different  SOCs  on  the  same  cell  (cell  1). 
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Fig.  9.  Voltage  estimation  results  of  ECM  and  EECM  under  the  parameterization  profile  Voltage  estimation  results  of  ECM  and  EECM  under  the  parameterization 

on  the  same  cell  (cell  1):  (a)  at  50%  SOC;  (b)  at  5%  SOC.  Profile  on  a  different  cell  (cell  2):  (a)  at  50%  SOC;  (b)  at  5%  SOC. 


4.  Model  verification  under  different  load  profiles 

In  this  section,  the  voltage  estimation  accuracy  of  the  EECM  is 
verified  under  different  load  profiles  based  on  the  identified  pa¬ 
rameters  in  Section  3.  To  validate  the  model  effectiveness  and  the 
advantage  over  the  original  ECM,  two  load  profiles  are 
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Fig.  10.  Voltage  RMSE  comparison  of  ECM  and  EECM  under  the  parameterization 
profile  at  different  SOCs  on  a  different  cell  (cell  2). 


implemented  in  the  verification  process,  namely  the  constant- 
current  discharge  pause  profile  (DPP)  in  Section  4.1  and  the  dy¬ 
namic  stress  test  (DST)  profile  in  Section  4.2.  It  is  to  be  noticed  that 
the  model  input  is  the  current  or  power  profile  in  the  whole  test 
period  indicating  an  offline  computing  process,  which  is  different 
from  the  online  SOC  estimation  case  in  Section  5  with  real-time 
voltage  update. 


4.1.  Results  under  constant- current  discharge  pause  profile 

The  performance  of  EECM  is  firstly  evaluated  under  a  set  of 
constant-current  discharge  pause  profiles  (DPP).  The  DPP  is  deter¬ 
mined  as  a  two-period  cycle,  in  which  the  battery  experiences 
constant-current  discharge  in  the  first  period  and  rests  in  the  sec¬ 
ond  period,  and  cycles  from  the  fully  charged  state  until  the 
discharge  voltage  limit.  The  long  constant-current  discharge  period 
aims  to  simulate  the  stable  driving  condition  with  relatively  con¬ 
stant  speed  and  power  demand,  and  the  rest  period  is  to  check  the 
voltage  estimation  accuracy  in  the  depolarization  process.  Three 
groups  of  DPP  tests  were  carried  out  with  the  current  values  10  A 
(0.5  C),  20  A  (1  C),  and  40  A  (2  C)  respectively,  in  which  the  lengths 
of  the  discharge  periods  were  respectively  20  min,  10  min,  and 
5  min  to  get  an  identical  discharge  coulomb  amount  and  hence  an 
identical  SOC  change  for  each  DPP  cycle,  in  order  to  examine  the 
influence  of  discharge  rate  on  model  accuracy  in  the  same  SOC  area. 
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Fig.  12.  Current  and  SOC  variation  of  discharge-pause  profile  (under  0.5  C  DPP). 


Each  test  was  with  a  uniform  20-min  rest  time.  The  current  profile 
and  the  corresponded  SOC  variation  are  illustrated  in  Fig.  12  with 
the  0.5C  DPP  as  an  example,  in  which  six  complete  DPP  cycles  are 
accomplished  and  the  battery  reaches  the  discharge  limit  in  the 
seventh  discharge  period.  It  is  to  be  noted  that  the  validation  ex¬ 
periments  in  Section  4.1  and  4.2  were  carried  out  on  cell  2  (with  the 
model  parameters  identified  on  cell  1)  to  check  the  model  appli¬ 
cability  considering  cell  inconsistency. 

Under  0.5  C  DPP  test,  the  estimated  voltages  of  the  ECM  and 
EECM  are  compared  with  the  measured  voltage  in  Fig.  13a,  with  the 
voltage  error  shown  in  Fig.  14a.  Note  that  the  error  is  determined  by 
the  measured  value  subtracted  by  the  estimated  value,  so  that  a 
positive  error  indicates  that  the  model  underestimates  the  voltage 
while  a  negative  error  implies  overestimation.  The  green  dashed 
lines  (in  the  web  version)  in  Fig.  14  reflect  the  10  mV  and  -10  mV 
error  boundaries.  It  is  shown  in  Figs.  13a  and  14a  that  both  the  ECM 
and  the  EECM  provide  accurate  results  at  the  beginning  and  middle 
periods  of  the  test  with  the  errors  within  ±10  mV.  Starting  from  the 
sixth  cycle  with  SOC  lower  than  19.5%,  the  estimated  voltages  show 
larger  deviation  from  the  real  value,  but  the  EECM  could  still  guar¬ 
antee  a  low  voltage  error  (lower  than  20  mV)  in  contrast  to  the  ECM 
with  maximum  60  mV  underestimation.  It  is  also  noticed  that  the 
voltage  accuracy  in  the  rest  period  is  improved  by  the  EECM,  which 
could  hence  be  implemented  in  the  adaptive  estimation  method  of 
open  circuit  voltage  [33  .  The  voltage  comparison  and  the  estima¬ 
tion  errors  under  1  C  DPP  test  are  demonstrated  in  Figs.  13b  and  14b, 
and  similarly  the  low-SOC-range  voltage  error  of  the  EECM  is  much 
lower  than  that  of  ECM  (max.  30  mV  versus  max.  100  mV).  The 
analysis  of  the  2  C  DPP  case  in  Figs.  13c  and  14c  also  leads  to  the  same 
conclusion  that  the  EECM  significantly  reduces  the  low-SOC-area 
voltage  error  (max.  50  mV  versus  max.  130  mV).  By  comparing  the 
voltage  error  under  different  discharge  currents  in  Fig.  14a,  b  and  c,  it 
is  clear  that  the  model  estimation  error  increases  at  the  current 
grows,  and  the  EECM  could  limit  the  voltage  error  in  a  much  smaller 
range  than  the  ECM  and  consequently  upgrade  the  low-SOC-area 
accuracy  under  different  discharge  currents. 

The  RMSE  values  under  different  SOC  ranges  and  currents  are 
then  compared,  in  which  the  voltage  RMSE  is  calculated  by  the 
mean  voltage  error  in  the  discharge  period  of  each  DPP  cycle.  Fig.  15 
lists  the  RMSE  values  of  EECM  and  ECM  related  to  SOC,  and  the  SOC 
value  is  treated  as  the  mean  SOC  of  each  complete  discharge  period, 
i.e.  92.0%,  75.9%,  59.8%,  43.7%,  27.6%,  and  11.5%  respectively.  The 
blue  markers  are  for  the  ECM  and  the  red  for  EECM,  while  the 
square,  circle,  and  triangle  markers  indicate  separately  the  errors  in 


Fig.  13.  Voltage  estimation  results  of  ECM  and  EECM  in  DPP  test  with  different  current 
values:  (a)  under  0.5  C  DPP;  (b)  under  1  C  DPP;  (c)  under  2  C  DPP. 


0.5  C,  1  C,  and  2  C  DPP  cases.  At  SOC  larger  than  20%,  the  RMSE 
values  are  relatively  independent  from  the  SOC  variation  for  both 
models.  The  lower-than-20%  SOC  leads  to  a  large  RMSE  increase, 
while  the  EECM-based  RMSE  values  are  reduced  by  more  than  50% 
compared  with  the  ECM  (e.g.  19  mV  versus  41  mV  in  2  C  DPP  case). 
Additionally,  the  ECM-based  RMSE  result  grows  rapidly  when  the 
current  increases,  while  the  EECM-based  error  is  well  controlled 
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Time  (min) 


Fig.  14.  Voltage  error  comparison  of  ECM  and  EECM  in  DPP  test:  (a)  under  0.5  C  DPP; 
(b)  under  1  C  DPP;  (c)  under  2  C  DPP. 


under  different  current  rates,  with  the  error  lower  than  10  mV  at 
high  SOC  and  under  20  mV  in  low-SOC  range. 

Different  from  the  traditional  ECM,  the  EECM  determines  the 
cell  potential  as  well  as  the  resistance  parameters  through  the 
surface  SOC  rather  than  the  macroscopic  SOC,  from  which  a  better 
reflection  of  the  nonlinearity  in  voltage  variation  is  resulted.  As  in 
Eq.  (26),  the  SOCsurf  must  be  determined  before  the  terminal 
voltage  calculation.  Taking  the  2  C  DPP  test  as  an  example,  the 
variation  of  the  macroscopic  SOC  and  the  SOCsurf  is  shown  in  Fig.  16, 


Time  (min) 


Fig.  16.  Variation  of  SOC  and  surface  SOC  (SOCsurf)  under  2  C  DPP. 


and  the  SOCsurf  is  always  smaller  than  the  average  SOC  due  to  a 
non-negative  discharge  current.  From  the  zoom  figure  of  the  last 
discharge  period,  it  is  illustrated  that  the  gap  between  SOC  and 
SOCsurf  increases  during  discharge  owing  to  the  progress  of  solid- 
phase  diffusion,  and  the  rest  minutes  diminish  the  SOC-SOCsurf 
gap  indicating  the  balancing  process  of  the  surface-bulk  concen¬ 
tration  difference  in  the  electrode  particle. 

Under  different  current  rates,  the  SOC-SOCSUrf  gap  in  relation 
to  SOC  is  shown  in  Fig.  17.  As  expressed  in  Fig.  4,  the  SOC-SOCsurf 
deviation  (ASOC)  endows  the  EECM  with  nonlinear  capability  in 
AUSurf  determination  instead  of  the  linear  voltage  approximation 
from  the  macroscopic  SOC  point,  which  results  in  the  advantage 
over  the  ECM  in  voltage  estimation  performance.  Under  larger 
discharge  current,  the  ASOC  is  larger  for  the  EECM  and  could 
consequently  represent  the  voltage  nonlinearity  more  evidently, 
and  the  model  benefit  over  the  ECM  is  for  that  reason  larger, 
which  is  in  accordance  with  the  RMSE  results  in  Fig.  15.  The  ASOC 
tends  to  increase  as  the  discharge  continues,  and  the  ASOC-vari- 
ation  is  probably  caused  by  the  variation  of  the  identified  DSOC- 
related  parameters  (kSD,i,  tSd,i,  /<sd,2,  and  tSd,2)  considering  SOC 
change.  To  sum  up,  the  performance  of  the  EECM  is  validated 
under  DPP  tests  with  different  rates,  and  the  advantage  over  the 


Fig.  17.  Variation  of  the  SOC-SOCsurf  deviation  (ASOC)  in  relation  to  SOC  in  DPP  test. 
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Fig.  18.  Power  profile  of  DST  test  cycle. 


ECM  is  shown  in  voltage  estimation  accuracy  especially  in  low- 
SOC  cases. 

4.2.  Results  under  DST  profile 

The  model  validation  is  then  carried  out  under  the  dynamic 
stress  test  (DST)  load  profile.  The  DST  profile  is  a  typical  dynamic 
experiment  cycle  to  test  the  battery  dynamic  properties  reflecting 
the  power  demand  in  dynamic  EV  driving  conditions  [34  .  The  DST 
is  here  realized  as  a  determined  power  profile,  with  the  largest 
discharge  power  200  W  to  keep  the  current  under  the  maximum 
allowable  discharge  current  80  A.  As  shown  in  Fig.  18,  an  experi¬ 
ment  cycle  is  determined  here  as  two  consecutive  DSTs  followed  by 
a  20-min  rest  period.  To  prevent  overcharging  at  the  fully-charged 
battery  state,  the  test  is  performed  starting  from  70%  SOC  until  the 
end  of  discharge. 

The  estimated  voltages  of  the  ECM  and  EECM  are  compared  with 
the  measured  voltage  under  DST  profile  in  Fig.  19.  A  total  of  ten 
complete  DST  cycles  are  accomplished,  and  the  4th  and  10th  cycles 
are  analyzed  in  detail.  Fig.  20a  corresponds  to  the  voltage  results  in 
the  4th  cycle  with  SOC  50.7%-43.9%  representing  the  middle-SOC 
range  and  Fig.  21a  illustrates  the  voltage  error  with  the  green 
dashed  lines  (in  the  web  version)  as  the  ±10  mV  boundaries.  Both 
the  ECM  and  the  EECM  provide  accurate  results  with  the  errors 


Fig.  19.  Voltage  estimation  results  of  ECM  and  EECM  under  DST  profile. 


Fig.  20.  Voltage  estimation  results  of  ECM  and  EECM  in  DST  test  cycles  with  different 
SOC  ranges:  (a)  4th  test  cycle  with  SOC  range  50.7%-43.9%;  (b)  10th  test  cycle  with  SOC 
range  9.7%-1.5%. 


Fig.  21.  Voltage  error  comparison  of  ECM  and  EECM  in  DST  test  cycles  with  different 
SOC  ranges:  (a)  4th  test  cycle  with  SOC  range  50.7%-43.9%;  (b)  10th  test  cycle  with  SOC 
range  9.7%-1.5%. 
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Fig.  22.  Voltage  RMSE  comparison  of  ECM  and  EECM  under  DST  profile  at  different 
SOCs. 

within  ±10  mV  and  fulfill  the  normal  onboard  accuracy  require¬ 
ment  [35]  while  the  EECM-based  error  is  slightly  smaller.  Figs.  20b 
and  21b  are  related  to  the  10th  cycle  in  low-SOC  area  with  SOC 
9.7%-1.5%,  and  the  EECM  presents  obviously  better  voltage  simu¬ 
lation  results  with  the  RMSE  decreased  by  64%  (17.9  mV  versus 
49.6  mV)  and  the  maximum  error  decreased  by  66%  (52  mV  versus 
151  mV).  As  a  result,  the  low-SOC-range  voltage  error  under  dy¬ 
namic  load  conditions  could  be  significantly  reduced  with  the 
implementation  of  the  EECM. 

The  RMSE  values  under  different  SOC  ranges  are  then  compared 
in  Fig.  22  with  the  voltage  RMSE  calculated  by  the  mean  error  in 
two  consecutive  DSTs  of  each  cycle,  in  which  the  SOC  value  is 
treated  as  the  mean  SOCs  of  each  consecutive  DST  period,  i.e.  66.6%, 
60.3%,  53.7%,  47.4%,  40.8%,  34.1%,  27.3%,  20.4%,  13.3%,  and  5.8% 
respectively.  At  SOC  larger  than  20%,  the  RMSE  values  are  relatively 
stable  and  both  the  ECM-based  and  EECM-based  errors  are  kept 
under  10  mV.  For  the  ECM,  the  RMSE  increases  largely  at  SOC  lower 
than  20%,  especially  in  the  lower-than-10%  SOC  area,  while  the 
EECM-based  error  is  well  controlled  under  20  mV.  It  is  verified  that 
although  the  EECM  could  only  slightly  increase  the  voltage  accu¬ 
racy  in  high-SOC  area,  the  effectiveness  of  the  extended  ECM  is 
explicit  for  the  error  reduction  in  low-SOC  cases. 


Fig.  23  shows  the  variation  of  the  macroscopic  SOC  and  SOCsurf 
in  the  EECM  under  DST  profile.  From  the  zoom  figure  of  the  last  test 
period,  it  is  known  that  the  gap  between  SOC  and  SOCsurf  (i.e.  ASOC) 
varies  frequently  due  to  the  dynamic  discharging  and  charging 
loads  in  DST.  In  the  high-power  discharge  periods  with  maximum 
current  72  A  (3.6  C),  the  ASOC  increases  rapidly  and  reaches  an  apex 
of  1.5%,  which  makes  the  voltage  nonlinearity  more  obviously  and 
results  in  a  larger  advantage  over  ECM,  as  reflected  by  the  error 
curves  in  Fig.  21b. 

5.  SOC  estimation  results  based  on  the  extended  equivalent 
circuit  model 

5.2.  SOC  estimation  system  based  on  the  EECM  through  extended 
Kalman  filtering 

In  the  vehicle  onboard  battery  management  system,  the  bat¬ 
tery  model  is  the  basis  of  state  estimation,  which  includes  the 
voltage  model,  the  aging  model,  the  temperature  model,  and  the 
safety  model,  etc.  to  determine  the  battery  remaining  discharge 
energy  (Erde),  the  power  capability,  the  remaining  useful  life,  and 
the  safety  state  [36-38  .  The  battery  voltage  model  discussed  in 
this  paper  is  mainly  employed  to  estimate  the  state  of  charge 
(SOC)  real-timely,  which  is  further  used  for  the  remaining 
discharge  energy  (Erde)  and  remaining  driving  range  (RDR) 
determination.  There  exist  a  large  number  of  model-based 
methods  for  online  SOC  estimation  13,36,39-41],  in  which  a  se¬ 
ries  of  adaptive  algorithms  (i.e.  Kalman  filtering,  extended  Kalman 
filtering,  sigma-point  Kalman  filtering,  particle  filtering,  H-in- 
finity,  etc.)  attract  high  attention.  These  methods  combine  the 
coulomb  counting  results  and  the  model  voltage  outputs  to  obtain 
good  accuracy.  The  extended  Kalman  filtering  (EKF)  is  applicable 
for  nonlinear  system  and  the  computational  process  is  relatively 
simple  compared  with  most  other  adaptive  methods,  so  that  it  can 
be  executed  very  accurately  and  continuously  in  onboard  con¬ 
trollers.  There  exists  a  considerable  influence  on  EKF  performance 
from  the  accuracy  of  the  employed  model,  and  the  model  accuracy 
could  hence  be  evaluated  through  the  SOC  estimation  results 
based  on  different  models.  As  a  result,  the  EKF  estimator  is 
respectively  carried  out  based  on  ECM  and  EECM  to  obtain  the  SOC 
estimation  outcomes. 

The  battery  model  needs  to  be  converted  into  a  discrete  state- 
space  form  in  real-time  applications.  At  the  present  time  point  t/<, 
the  state  equation  and  the  measurement  equation  of  the  battery 
system  are  expressed  in  Eqs.  (28)  and  (29),  in  which  x /<  is  the  system 
state  vector.  The  terms  Uk  and  yj<  represent  respectively  the  current 
input  and  output  which  stand  here  for  the  measured  current  /  and 
terminal  voltage  Uh  while  the  process  noise  and  the  measurement 
noise  are  written  as  w/<  and  v/<.  fixk-i,  Uk)  is  the  nonlinear  state 
transition  function  and  g(x/0  u/J  is  the  nonlinear  measurement 
function.  In  EKF,  both/(x/<_i,  u/<)  and  g(x/0  u/<)  are  linearized  in  each 
time  step  by  first-order  Taylor  expansion  and  hence  expressed  as 
the  parameter  matrices  Aj0  %  Q,  and  D/<,  as  in  Eq.  (30).  Since  the 
battery  parameters  are  affected  by  state  factors  such  as  SOC  and 
temperature,  the  parameters  in  the  matrices  need  to  be  looked  up 
during  operation,  indicating  a  linear-parameter-varying  (LPV)  sys¬ 
tem.  The  coulomb-counting  information  combined  in  the  EKF 
estimator  is  shown  in  Eq.  (31 ),  which  is  used  to  determine  the  prior 
state  estimate  x ^  in  a  new  time  step.  u/<  is  the  current  input  at  time 
k,  At  is  the  sample  time,  Qst  is  the  battery  standard  capacity,  and  ri  is 
the  coulomb  efficiency  which  is  determined  as  1  for  the  later 
calculation. 

*k  =/(Xk-l,Ufc)  +  Wk 


(28) 
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Yk  =  g(*io  u;<)  +  vk 


(29) 


Ai 


G 


d*k- 1 

Sg(xk,uk) 


x=  Xk_! 


-  Bk 


dxk 


-  ,  Dk 
x=  xk 


duk 

dg(Xk ,  Ufc) 


X=  x/c-l 


(30) 


x=  Xk 


^ECM,/<  — diag(l  ,0(r; po,l  )  >#(r po,2)) 

—At 


^ECM,k~  (ostr;-3600’^p0'1  (”*  ^(TP°.l(J’^P°.2^  ^(rpo,2))^ 


CECM,fc=  (^(SOCk)  -1,-1 


^ECM,k— _^o 


(33) 


SOC/f  SOCfc_!  ukQsvV.360Q  (31) 

The  SOC  estimator  based  on  the  two-state  ECM  is  firstly 
expressed  with  the  state  vector  xEcm,/<  in  Eq-  (32),  in  which  SOQ, 
Upoxk,  and  Up0t2,k  are  the  SOC  value  and  the  polarization  voltages 
on  two  RC  components,  respectively.  The  parameter  matrices 
AecmJo  £>ecm jo  Cecm jo  and  Decm,/<  are  shown  in  Eq.  based  on 
Eqs.  (18)  and  (19),  in  which  the  term  #(  •)  is  the  simplification  of 
the  exponential  relationship  as  in  Eq.  (34).  Based  on  the  identi¬ 
fication  results  in  Section  3,  the  parameters  in  the  matrices  are 
looked  up  in  relation  to  SOC,  in  which  the  term  D0cv/S OC(SOC/J 
in  CEcm,/<  is  the  ratio  of  Docv  and  SOC  at  the  present  macroscopic 
SOC. 


6(ri)  =  exp(  -  A t/Tj)  (34) 

Then  the  EECM-based  SOC  estimator  is  determined  with  the 
state  vector  xEEcm,/<  in  Eq.  (35),  in  which  SOCsurfi,<,  UDLXk,  Udl,2 jo 
AS0Ci  /o  and  ASO C2j<  are  respectively  the  surface  SOC  value,  the 
polarization  voltages  of  the  double-layer  process,  and  the  SOCSUrf 
offsets  due  to  solid-phase  diffusion.  The  macroscopic  SOC  is 
accessible  from  xeecm,/<  in  Eq.  (36).  Eq.  (37)  shows  the  parameter 
matrices  Aeecmjo  Beecmm,  Ceecm jo  and  Deecm,/<  based  on  Eqs. 
(23)-(27  .  Different  from  the  ECM,  the  varying  parameters  in  the 
EECM  matrices  are  looked  up  by  SOCsurf  value,  and  the  term  D0cv/ 
SOC(SOCSUrf,/<)  in  CEecm,i<  is  the  ratio  of  D0cv(SOCSUrf,/<)  and  the  pre¬ 
sent  surface  SOC  SOCsurf;/<. 

xEECM,/<  =  (SOC^*,  ^DL,l,/o  ^DL,2,/o  ^SOC} *,  ASOC2  /^  (35) 


^ECM,/<  —  (S0C/o  Dp0  1  /0  Up0  2Ji 


(32)  SOC/c—  SOCsurf  x  +  ASOCv<  +  ASOC2  * 


(36) 


/I 


A 


EECM,/< 


l-#(rSD,l)  ^-^(TSD,2)\ 
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#(TDL,2) 
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^(rSD,2)  / 
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Fig.  24.  Schematic  of  SOC  estimation  process  through  extended  Kalman  filter. 


The  iteration  process  of  the  EKF  SOC  estimator  is  shown  in 
Fig.  24  based  on  the  determined  state  spaces.  According  to  the  state 
xk_^  and  the  error  covariance  P+_t  of  the  previous  time  step  as  well 
as  the  input  ii/<  of  the  present  time  point,  the  SOC  estimator  derives 
the  present  post  state  estimate  xk  and  the  error  covariance  esti¬ 
mate  P^,  in  which  Q/<  and  P/<  are  respectively  the  process  and 
measurement  noise  covariances  and  I<k  is  the  calculated  Kalman 
gain.  In  the  iteration  process,  the  prior  state  xk  and  the  corre¬ 
sponded  model  voltage  output  yk  are  calculated  at  first,  and  the 
error  e /<  between  the  measured  and  model-estimated  voltages  is 
then  combined  with  the  Kalman  gain  I< f  to  update  the  state  vector. 

5.2.  Results  and  comparison 

The  SOC  results  of  the  EKF  estimator  based  on  ECM  and  EECM 
are  compared  in  this  section,  in  which  the  experiment  data  under 
DST  profile  is  implemented  and  the  sample  time  At  is  set  as  0.1  s. 
For  the  initial  state  xo,  the  initial  SOC  (or  SOCsurf  for  the  EECM)  is 
determined  as  100%,  i.e.  30%  larger  than  the  real  value  70%,  and  the 


Fig.  25.  SOC  estimation  results  comparison  through  EKF  SOC  estimator  based  on  ECM 
and  EECM  in  DST  test. 


other  voltage  terms  and  SOCsurf  offsets  in  the  state  vector  (Lfpo,i,o 
and  [7po,2,o  for  ECM,  and  L1Dl,i,o,  ^dl,2,o.  ASOCii0,  and  ASOC2,0  for 
EECM)  are  defined  as  0.  It  is  noticed  that  the  SOC  updating  result  in 
EKF  is  affected  by  the  initial  error  covariance  Po  as  well  as  the  noise 
covariances  Q/<  and  Rj0  and  the  corresponded  terms  in  Po,  Qj<}  and  P/< 
of  the  ECM-based  and  EECM-based  estimation  systems  are  hence 
determined  as  identical  to  avoid  the  influence  of  noise  difference. 

Fig.  25  illustrates  the  SOC  estimation  results  based  on  the  two 
models,  in  which  the  reference  SOC  is  determined  by  coulomb 
integration  with  an  accurate  initial  state.  It  is  shown  that  both 
models  converge  rapidly  towards  the  real  SOC  to  correct  the  initial 


Fig.  26.  SOC  error  comparison  of  ECM  and  EECM  at  different  SOC  ranges:  (a)  complete 
test  process:  (b)  SOC  range  9.7%-1.5%. 
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error,  and  the  SOC  errors  in  a  wide  SOC  range  are  relatively  small. 
The  error  increase  of  ECM-based  SOC  is  noticed  in  the  last  test 
period  with  SOC  9.7%— 1.5%,  and  the  error  grows  rapidly  especially 
under  high  discharge  rate  corresponding  to  the  SOC-rapid- 
decreasing  period  in  the  figure.  In  contrast,  the  EECM-based  SOC 
maintains  its  accuracy  in  low-SOC  area  and  the  performance  is  little 
influenced  by  the  discharge  rate. 

The  error  between  the  model-estimated  SOC  and  the  real  SOC  is 
shown  in  Fig.  26,  in  which  the  error  is  determined  by  the  real  SOC 
subtracted  by  the  estimated  value,  and  a  positive  error  indicates  an 
underestimation  of  SOC  while  a  negative  error  implies  SOC  over¬ 
estimation.  Both  models  show  good  accuracy  in  a  wide  SOC  area 
with  the  SOC  error  within  1%,  and  the  ECM-based  result  becomes 
worse  at  low  SOC  with  a  maximum  error  of  2%  while  the  EECM 
always  presents  very  limited  errors  under  0.5%.  Consequently,  the 
EECM  helps  to  improve  the  estimation  accuracy  in  low-SOC  range 
evidently. 

Known  from  the  analysis  in  Section  4.2,  the  ECM-based  voltage 
error  reaches  a  maximum  of  150  mV  in  low-SOC  area,  while  the 
corresponded  largest  SOC  error  is  2%  and  seems  not  so  serious.  It  is 
probably  because  that  the  SOC-Uocv  curve  is  steeper  at  low  SOC 
and  hence  leads  to  a  relatively  smaller  SOC  estimation  error  with 
respect  to  the  same  voltage  error.  It  seems  that  the  ECM-based  SOC 
error  would  not  influence  the  Erde  calculation  largely,  but  since  the 
battery  model  is  also  used  to  predict  the  future  terminal  voltage  in 
the  Erde  determination  process  as  in  Eq.  (1 ),  the  large  voltage  error 
together  with  the  SOC  error  would  lead  to  larger  mistakes  in  Erde 
prediction  results.  The  EECM  is  comparatively  with  a  large  advan¬ 
tage  in  the  ERD e  and  RDR  determination. 

Additionally,  the  SOC  error  based  on  ECM  tends  to  be  negative 
reflecting  that  the  ECM  tends  to  overestimate  the  SOC  during 
operation,  while  a  positive  error  is  normally  resulted  for  the  EECM 
indicating  underestimation.  In  the  PEV  case,  the  SOC  over¬ 
estimation  implies  that  the  battery  energy  would  be  exhausted 
earlier  than  expected  and  the  vehicle  would  operate  in  the  limp- 
home  mode  before  arriving  at  the  destination,  which  increases 
the  range  anxiety  of  PEV  passengers  [2  .  On  the  contrary,  the  SOC 
underestimation  for  the  EECM  indicates  more  battery  energy  and 
more  remaining  range  than  expectation,  which  is  relatively  safe 
and  could  help  reduce  the  users'  range  anxiety  for  electric  vehicles. 

Fig.  27  compares  the  SOC  RMSE  under  different  SOC  ranges  for 
the  ECM-based  and  EECM-based  EKF  estimators,  and  the  RMSE  is 
determined  in  the  periods  identical  as  in  Section  4.2.  At  SOC  larger 
than  20%,  the  estimation  accuracy  of  both  the  ECM  and  the  EECM  is 


Fig.  27.  SOC  RMSE  comparison  of  ECM  and  EECM  at  different  SOCs. 


satisfying  with  errors  under  0.4%.  For  the  SOC  under  20%  and 
especially  for  the  lower-than-10%  case,  the  RMSE  for  the  ECM  in¬ 
creases  largely  to  a  maximum  of  1%,  while  the  EECM-based  error  is 
still  well  regulated  under  0.4%.  As  a  result,  the  EECM  could 
evidently  enhance  the  SOC  estimation  accuracy  in  low-SOC  range. 

6.  Conclusion 

In  this  work,  an  extended  equivalent  circuit  model  (EECM)  is 
presented  with  enhanced  performance  in  low  state-of-charge 
(SOC)  area  compared  with  the  traditional  equivalent  circuit 
model  (ECM).  The  EECM  represents  the  battery  electrochemical 
process  with  electrical  components,  in  which  two  RC  components 
reflect  the  double  layer  effect  and  two  one-state  elements  (DSOCs) 
indicate  the  solid-phase  diffusion.  Different  from  the  ECM  which 
reflects  the  solid-phase  diffusion  by  polarization  voltage,  the  two 
DSOC  elements  in  the  EECM  calculate  the  difference  between  the 
macroscopic  SOC  and  the  particle  surface  SOC  (SOCsurf),  and  the 
obtained  SOCsurf  is  then  implemented  to  determine  the  terminal 
voltage.  The  EECM  exhibits  better  nonlinear  performance  thanks  to 
the  SOCsurf  introduction  and  could  hence  provide  satisfactory 
estimation  accuracy  in  low-SOC  range. 

The  effectiveness  of  the  EECM  is  validated  through  comparison 
with  the  commonly-used  two-state  ECM  on  a  type  of  prismatic 
NCM-G  battery.  The  genetic  algorithm  is  employed  in  model 
parameter  identification,  in  which  the  EECM  presents  a  better 
fitting  result  than  the  ECM  at  SOC  lower  than  20%.  The  voltage 
estimation  accuracy  of  ECM  and  EECM  is  then  compared  under  two 
types  of  load  profiles  in  SOC  range  0-100%.  Under  three  groups  of 
constant-current  discharge  pulse  profiles  (DPP)  with  0.5  C,  1  C,  and 
2  C  current  rates,  the  voltage  root  mean  square  errors  (RMSEs)  in 
low-SOC  area  are  all  reduced  by  more  than  50%  through  the  EECM 
implementation.  The  error  increase  under  larger  current  is  also 
smaller  for  the  EECM,  indicating  its  applicability  for  the  high-rate 
cases.  Under  DST  profiles,  the  voltage  RMSE  in  low-SOC  case  is 
reduced  by  64%  based  on  the  extended  model,  implying  the  effec¬ 
tiveness  in  dynamic  operating  conditions.  For  SOC  estimation,  the 
performance  of  EECM  and  ECM  is  discussed  through  the  extended 
Kalman  filtering  (EKF)  estimator.  The  result  based  on  the  EECM 
shows  a  significant  advantage  in  low-SOC  range  over  the  traditional 
ECM,  while  the  EECM-based  SOC  error  is  limited  under  0.5%  in  the 
whole  SOC  range.  Since  the  model  complexity  is  well  controlled, 
the  EECM  could  be  onboard  implemented  without  difficulty  to 
obtain  a  preciser  SOC  and  voltage  estimation  result  in  the  entire 
SOC  range,  so  as  to  improve  the  prediction  accuracy  of  the 
remaining  battery  energy  and  the  remaining  driving  range  in  EV 
applications.  The  future  work  will  focus  on  the  possible  simplifi¬ 
cation  of  the  EECM  and  the  model  performance  validation  under 
different  temperatures  and  aging  conditions. 
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Abbreviations,  nomenclature,  and  subscripts 

Abbreviations 
Ah  ampere  hour 

BMS  battery  management  system 

CPE  constant  phase  element 
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DPP  discharge  pause  profile 

DST  dynamic  stress  test 

ECM  equivalent  circuit  model 

EECM  extended  equivalent  circuit  model 

EIS  electrochemical  impedance  spectroscopy 

EKF  extended  Kalman  filtering 

GA  genetic  algorithm 

HPPC  hybrid  pulse  power  characterization 

LPV  linear  parameter  varying 

NCM  LiNixCoyMni_x_y02 

PEV  pure  electric  vehicle 

P2DM  pseudo  two  dimensional  model 

RDR  remaining  driving  range 

RMSE  root  mean  square  error 

SEI  solid  electrolyte  interface 

SOC  state  of  charge 

SPM  single  particle  model 

Nomenclature 

A  electrode  plate  area  (m2) 

Cb  bulk  average  lithium  concentration  of  electrode  particle 

(mol  nrr3) 

csurf  lithium  concentration  at  electrode  particle  surface 
(mol  m~3) 

Cb,o%  bulk  average  lithium  concentration  of  electrode  particle 
at  SOC  =  0  (mol  nrT3) 

Ds  solid  phase  diffusion  coefficient  (m2  s-1) 

Erde  battery  remaining  discharge  energy  (Wh) 

j  current  density  (A  rrT3) 

j'sd  solid-phase-diffusional  current  density  (A  m~3) 

/<sd,c  concentration  difference  parameter  in  solid  diffusion 
(mol  itT3  A-1) 

/<sd  SOC-difference  parameters  in  solid  diffusion  (m3  mor1) 
Qnom  nominal  discharge  capacity  (Ah) 

Qst  standard  discharge  capacity  (Ah) 

R0  Ohmic  resistance  (Q) 

Rs  radius  of  spherical  particles  (m) 

SOCsurf  surface  state  of  charge  (-) 

ASOC  difference  of  SOC  and  SOCsurf  (-) 

Docv  open  circuit  voltage  (V) 

D0  instantaneous  voltage  drop  in  discharge  (V) 

Dsei  voltage  drop  on  SEI  layer  (V) 

Ut  terminal  voltage  (V) 

5  electrode  or  separator  thickness  (m) 

r  time  constant  (s) 

£s  active  material  volume  fraction 

0(t)  exponential  relationship:  exp (-t/r) 

Subscripts 

Avg  average  value 

DL  related  to  double  layer 

Lim  related  to  values  at  discharge  voltage  limit 

N  related  to  negative  electrode 


P  related  to  positive  electrode 

po  related  to  polarization  terms  in  ECM 
SD  related  to  solid-phase  diffusion 

surf  related  to  values  at  particle  surface 
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